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Abstract. Time shifting the outputs of Gravitational Wave detectors operating in 
coincidence is a convenient way to estimate the background in a search for short 
duration signals. However this procedure is limited as increasing indefinitely the 
number of time shifts does not provide better estimates. We show that the false 
alarm rate estimation error saturates with the number of time shifts. In particular, for 
detectors with very different trigger rates this error saturates at a large value. Explicit 
computations are done for 2 detectors, and for 3 detectors where the detection statistic 
relies on the logical "OR" of the coincidences of the 3 couples in the network. 
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1. Introduction 

Kilometric interferometric Gravitational Wave (GW) detectors such as LIGO [TJ or 
Virgo [2] have been taking data with increasing sensitivities over the past years 
[4], [5J, [6l [Tj, [HI E] - It is expected that short duration GW events, e.g. the so-called bursts 
emitted by gravitational collapses or the signals emitted by compact binary inspirals, 
are very rare. Moreover the output of the detectors is primarily (non Gaussian) noise, 
and this background noise is in general not modeled. This implies that with a single GW 
detector it is very difficult to estimate the background event rate, and then to assess 
the significance of some GW candidate. 

On the contrary when dealing with a network of detectors (that means in practice 
at least two detectors of the same class), there is a conventional and simple way to 
estimate the background, that is the rate of coincident events due to detector noise. 
This consists of time shifting the search algorithm outputs (or triggers) of each detector 
with respect to the other (s), by some unphysical delays, much larger than the light travel 
time between the detectors and much larger than the typical duration of an expected 
GW signal. The next step is then to look for coincidences between shifted triggers just 
as if the shifted streams were synchronized. As we deal with a priori rare events, we 
need to set in practice low false alarm rates in the analysis. The question then arises of 
how many time slides are needed for correctly estimating the background and especially 
its tails where rare (non-Gaussian noise) events lay. Note, that in practice in burst 
or binary inspiral searches, a hundred or more time slides are done [21 El [7], due in 
particular to limited computational resources. Such a limitation of course depends on 
the duration of the different detectors data streams and on the complexity of consistency 
tests performed on coincident triggers. For example, time slides computation for one 
year of data sampled at 16384 Hz (LIGO) or 20 kHz (Virgo) can rapidly become a 
computational burden, especially when computationally intensive consistency tests like 
the x 2 ve to [10] are used. 

In this paper, we show that the precision on the background estimation, using time 
slides of trigger streams, is in fact limited. Indeed the variance of the false alarm rate 
estimation does not indefinitely decrease as the number of time slides increases as we 
would naively believe. On the contrary this variance saturates at some point, depending 
on the trigger rates chosen in each individual detector and on the coincidence time 
window set for identifying coincident events in the network of detectors. 

After introducing the general definitions in section [2j we give explicit formulas for 
the two-detector and the three-detector case in respectively Sec. [3] and [U In the latter 
we restrict ourselves to the particular analysis scheme where we are looking at the union 
(logical "OR") of the three couples of detectors. This is actually the configuration of 
interest, since it is more sensitive than simply searching for triple coincidences (logical 
"AND") [IT] . In each case (2 or 3 detectors) we check the analytical result with a Monte 
Carlo simulation and find excellent agreement. In section [5] these results are applied 
and discussed using typical parameters of a GW data analysis. 



Background estimation by time slides 



3 



2. Definitions 



2.1. Poisson approximation for trigger generation 

Background triggers in the detectors are due to rare glitches. Often these glitches come 
in groups, but most analysis pipelines cluster their triggers, so each glitch group results 
in only one final trigger. This clustering procedure is reasonable as long as the resulting 
trigger rate is much lower than the inverse of the typical clustering time length. In this 
limit the clustered triggers are independent events. Thus, throughout this paper we 
will assume that each detector produces random background triggers, which are Poisson 
distributed in time. 



2.2. Problem description 

We look then at the coincidence between two Poisson processes. The single 
interferometer trigger rate will be noted FA 1? FA 2 , the coincidence rate will be 
simply noted FA. We denote by FA(T) the rate resulting from counting the number of 
coincidence between two data streams, that are shifted by some time T. In particular 
for zero lag (T = 0) the measured rate is FA(0). So the quantities with tildes are the 
experimentally measured rates, and the quantities without tildes are the actual Poisson 
distribution parameters. The purpose of the paper is to study the properties of the time 
shifting method, which uses FA = J2k=i FA(7fc) as an estimator of FA, where R is the 
number of time slides and 7^ is the k th time slide stride. 



2.3. Poisson process model 

To model a Poisson process with event rate FAi we discretize the data stream duration 
T with bins of length At, the discretization time scale, e.g. either the detector sampling 
rate or the clustering time scale. Thus, for each bin an event is present with a probability 
p = FAi At Q. 

To ease the calculation we describe the Poisson process realizations with a 
continuous random variable. We take x uniformly distributed in the volume [0, 1]^ 
where N = £ is the number of samples, then compare Xk (the k th coordinate of x) 
with p. When x^, < p there is an event in time bin k, otherwise there is none. Thus 
x characterizes one realization of a Poisson process, and it can be easily seen that the 
uniform distribution of x leads to a Poisson distribution of events. 



2.4- Coincidences 

We choose the sampling At to be equal to twice the coincidence time window r c , in order 
to simplify the modeling of the coincidence between two processes. More precisely, for 
two Poisson processes with event rates respectively FAi and FA 2 , we define a coincidence 

| Here we model the Poisson process by a binomial distribution, recalling that when p <C 1 the binomial 
distribution tends toward a Poisson distribution 
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when there is an event in the same time bin k for both processes. This is different from 
the usual definition, where events are said in coincidence when they are less than a time 
window apart. This binning time coincidence has on average the same effects as defining 
as coincident events that are less than ±|Ai = ±r c apart. The analytical results are 
derived using this non standard definition, but they are in precise agreement with Monte 
Carlo simulations that are performed using the usual definition of time coincidence. 



3. The case of two detectors 

3.1. Time slides between two detectors 

Let x, y G [0, 1]^ be two realizations of Poisson processes with respectively p = FA 1 At 
and q = FA 2 At. There is a coincident event in time bin k when x k < p and y k < q. So 
the total number of coincidences for this realization is 

N 

^l(x k <p)l(y k <q) (1) 
k=i 

where 

1(a) = 1 if a is true 

1(a) = if a is false 

Thus the mean number of coincidences without time slides is as expected 

N 

Mean =/•••/ /••• / ^ Xk < p) 1 ^ < <?) dzi • • • dx N dyi ■ ■ ■ dyjy = Npq. 

J xi Jx N Jyi Jy N k=1 ^ ' 

(3) 

To consider a number R of time slides we take a set of R circular permutations of 
[1, N\. Time-sliding a vector x by the circular permutation n transforms the vector x 
into the vector of coordinates x n ^)- Then the mean number of coincidences is simply 



Mean 



/••• / /••• / ^ < < ^ dV = N py> ( 4 ) 

J XI Jxn Jw Jvn _ i. 



'X! Jx N Jyx Jy N n k 

thus there is no bias resulting from the use of time slides. 



3.2. Computation of the variance 



In order to have an estimate of the statistical error, we compute the variance with R 
time slides. The second moment is 



M, 



/'■■/ /"'/ ^J2J2^ Xk < P) 1 (Vn(k) < q) 

Jxi Jx N Jyi Jy N - ft n k 



7T1 7T2 k 




Figure 1. The solid line is the analytical formula ([7]) of the variance and the dashed 
line is the Monte Carlo variance as a function of the number of times slides. The Monte 
Carlo has been performed with FAj = 0.7 Hz, FA2 = 0.8 Hz, r c = 1ms, 500 trials and 
aT — 10 s data stream length. 



We can then exchange integrals and sums. To compute the integrals we distinguish 
two cases: when k ^ I the integrals on x^ and x\ are independent, and the integration 
over Xi, . . . ,xjy gives a p 2 contribution; otherwise the integration gives a p contribution. 
Analogously for the y variables we obtain q or q 2 depending on whether 7Ti(/c) = ^(Z) 
or not. 



The computation of this integral, detailed in Appendix A, yields 



Var 



Npq 
Npq 



R +P 



q + 



pq-(p + q) 
R 



2pq 



R +P + q 



(6) 
(7) 



where the last line is an approximation in the limit p, q, -5 <C 1, which is reasonable as 
far as GW analysis is concerned. 



3.3. Interpretation 

Each term in equation Q can be interpreted. The is what we would expect if 
we considered R independent Poisson process realizations instead of R time slides. The 
p+q comes from the estimation of the Poisson process event rate. Indeed, the estimation 
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of the event probability p from a single realization of a Poisson process with a mean 
number of events Np isp = p + 5p, where 5p is the random statistical error with variance 
(5p 2 ) = jj. This yields the mean rate of coincidences 

Mean = Npq = N (p + dp) (q + 8q) ~ Npq + Npdq + Nqdp (8) 

which corresponds to a variance of {N 2 p 2 Sq 2 + Nq 2 5p 2 ) = Npq{p + q), because 5p and 
5q are independent errors. Thus, when using only one realization for the single detector 
triggers, we have a statistical error on the single detector process rate. This statistical 
error is systematically propagated to the coincidence rate of each time slide, that yields 
the extra terms in the variance as compared to independent process realizations. One can 
see that this extra term is important when -| < max(jo, q); for cases where the coincident 
false alarm rate is maintained fixed (pq constant), the effect is most noticeable when p 
and q are very different. 

This gives an estimate of the variance of the number of coincident events in a data 
stream of length T. After converting to the estimation of the coincidence false alarm 
rate we obtain 



Mean^ = = FAiFA 2 At, (9) 



Mean 

¥ 

Var . . At 

Var a = Jv~ FAiFA 2 — 



4 + FAiAt + FA 2 At 
K 



(10) 



To verify these results, a Monte Carlo simulation has been performed. The Poisson 
processes are created as described in section El using a sampling rate of 16384 Hz, then a 
simple coincidence test with a window of r c = 1 ms is applied. The time shifts are done 
by adding an integer number of seconds to all events and applying a modulo T operation. 
The formula has been tested using 500 realizations of T = 10 4 s long Poisson processes, 
and using between 1 and 10 4 time slides for each realization. Figure Q] shows that the 
analytical formula ([3) and the Monte Carlo agree well for any number of time slides, and 
that the variance starts saturating when a few hundred time slides are used. We can 
see that the identification of the sampling time and the coincidence time window has no 
consequence on the result, the choice between binning and windowing the coincidences 
is a higher order effect. 



3.4- Straightforward extensions of the model 

In real data analysis, there are times when one of the detectors does not take science 
quality data for technical reasons. Thus, the data set is divided into disjoint segments, 
and the background estimation is often done by circular time slides on each segment 
separately. Afterwards the results from all the segments are combined to get the 
background false alarm estimation. The computation discussed above extends to this 
case with minimal changes. The circular permutations have to be changed to circular 
by block permutations, everything else being kept identical. 



Background estimation by time slides 



7 



Another caveat is that for real data analysis the coincidence procedure is often 
more complicated. Some of those complications are event consistency tests, e.g. do the 
two coincident events have a similar frequency? We can model this by adding some 
parameter / distributed uniformly in [0, 1] attached to each event, and then requesting 
a coincidence in the parameter /. 

For this model the results will be the same as those above, up to a factor of order 
1. Indeed, instead of applying a window of size At to our events, we are now working 
in a 2 dimensional (for instance time-frequency) space and using a rectangular window 
in this 2D parameter space. The procedure in both cases is the same — applying D 
dimensional rectangular windows to events distributed uniformly in a D dimensional 
space — up to the dimension of the space. 

4. The case of three detectors 

4-1. Time slides between three detectors 

In the case of three detectors one natural extension is to ask for events that are seen 
by at least two detectors, which means look for coincidence between two detectors for 
each detector pair, but counting the coincidences between three detectors only once. 
This "OR" strategy in a interferometer network has been shown to be more efficient 
than a direct three fold coincidence strategy ("AND" strategy) [TT]. For time slides, 
when shifting the events of the second detector with some permutation it, we also shift 
the events of the third detector by the same amount but in the opposite direction with 
tt" 1 . To write compact equations we abbreviate X = l(xk < p), Y = t(y n (k) < q), 
Z = l(^ 7r -i( fe ) < r), dV = dx\ ■ ■ -dx^dyi ■ ■ -dy^dzi ■ ■ -dzjy, where r = FA 3 At is the 
event probability per bin of the third detector and the vector z describes its realizations. 
Thus, the mean number of coincidences in the framework described in section [3Tl is 



Mean = / • • ■ / — ^ Yl ^ XY + YZ + XZ - 2XYZ] dV = N [pq + pr + qr - 2pqr] . 




k 



7T 



(11) 



4-2. Computation of the variance 



The second moment can be written compactly as 




[XYX'Y' + XZX'Z' + YZY'Z' + AXYZX'Y'Z' + 2XYX'Z' + 2XYY' Z' 

+2XZY'Z' - AXYX'Y'Z' - AXZX'Y'Z 1 - AYZX'Y'Z'} dV, (12) 



where the ' denotes whether the hidden variables are 7Ti, k or 7T2, I. 
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10 10 
number of time slides 



Figure 2. The solid line is the analytical formula l|T4|) of the variance and the dotted 
line is the Monte Carlo variance as a function of the number of time slides. The Monte 
Carlo has been performed with FAi = 0.04Hz, FA 2 = 0.08 Hz, FA 3 = 0.16 Hz, r c = 
31 ms, 500 trials and a T — 10 s data stream length. 



The computation of this integral, detailed in |Appendix B| yields 
N 



M 2 

R 



■ | (pq + pr + qr — 2pqr) 
+(R— 1) [pq(p + q + pq) + pr(p + r + pr) + qr(q + r + qr) + 6pqr — 4pqr(p + q + r)} 

+ [(R - 1)(N - 3) + (N - 1)] (pq + pr + qr - 2pgr) 2 j, (13) 

and can be approximated in the limit p, q, r, 1 by 

— + p + q + r H ). (14) 

ft pq + pr + qr J 

4-3. Interpretation 

Similarly to section 13.31 the extra terms in equation ( [141 1 can be explained through the 
error in the estimation of the single detector event rate. Using the same notations as in 
section 13.31 the mean coincidence number is 



Mean = iV (pq + pf + qr) ~ iV [pq + pr + qr + (q + r)5p + (p + r)5q + (p + q)5r] . 

(15) 
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Using the independence of estimation errors and recalling that (Sp 2 ) = -fe we obtain the 
variance of this mean value 

Var = N 2 [(5p 2 ){q + r) 2 + (5q 2 ){p + r) 2 + (5r 2 ){p + q) 2 ] 

— N [(pq + pr + qr){p + q + r) + 3pqr] , (16) 



that corresponds to the extra terms in equation (1T4I ). 

After converting to the estimation of the false alarm rate we obtain 



Mean^ ~ (FAiFA 2 + FA1FA3 + FA 2 FA 3 )Ai, (17) 

At 
~T 



Var^ ~ (FAiFA 2 + FA1FA3 + FA 2 FA 3 j 



I + FA lAf + FA.A, + FA 3 A f + FAiFAz 3 + F ™ 3 FA2FA3 *) - US) 

To check the 3 detector results we performed a Monte Carlo similar to the one of 
the 2 detector case (see section 13731 ) . The only difference is the number of detectors, and 
we choose a different coincidence window: r c = 31ms[§|. To check that the assumption 
of equal and opposite time slides does not influence the result, in the Monte Carlo the 
data in the second detector are shifted by and in the third detector by 37^. Figure 
[2] shows that the Monte Carlo and the 3 detector "OR" formula (PT4l) agree really well. 

4-4- The case of D detectors 

For the sake of completeness we can generalize the interpretation done in section 13.31 to 
the case of D detectors in the "AND" configuration. This generalization of equation (JSj) 
to D detectors yields a variance on the number of coincidences 

Var ~ NjlpJ i + I[Pj) » (19) 

i=l \ i=l jjti / 

where pi is the probability for detector i to have an event in a given time bin. 

The interpretation can also be generalized in the "OR" case, that is coincidence 
between any pair of detectors, although the computation is more cumbersome as detailed 
in 



Appendix C| and yields 



R ^ 2 

1=1 / ij^j,jj^k,kj^i 



(20) 



where pi is the event probability per bin in the i th detector. 

§ This accounts for the largest light travel time in the LIGO- Virgo network (27 ms) and some timing 
error in each detector. 
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5. Discussion 

We finally discuss the consequences of the above results on GW data analysis. To be 
able to put numbers into the equations we will look at a fiducial GW data taking run. 
We choose the run properties to be: 

• a duration of T = 10 7 s, that is roughly 4 months 

• two detectors with a light travel time separation of 25 ms, and we use the same time 
as the coincidence window, so that At = 50 mdjj], assuming perfect timing accuracy 
of trigger generators. 

• a desired coincidence false alarm rate of 10~ 8 Hz, i.e. one event every three years. 

We will look at two special cases of single detector threshold choice. One symmetric 
case, where thresholds are set so that the single detector trigger rate in each detector is 
roughly the same. One asymmetric case, where in one of the detectors there is only one 
trigger. This asymmetric case is extreme but instructive, because tuning the thresholds 
to obtain the best sensitivity often yields asymmetric trigger rates between different 
detectors. 

Symmetric detector case In this case we have the single detector trigger rate FAi = 
FA 2 = FA S = — 4.5 x 10~ 4 Hz, which gives using equation (fTOl) the fractional 
error of the false alarm estimation 

i 

^ ~ 3.2 ( — + 4.5 x 10~ 5 ) p = q = 2.25 x 10~ 5 

FA \R J 

for R = 100 
for R — > oo. 

So for 100 time slides we get a typical error of 30% in the false alarm estimation, 
and the error saturates at 2% for R > 20000. 

Asymmetric detector case In this extreme case the single detector trigger rates are 
FAi = ^ = 10~ 7 Hz and FA 2 = F ^ At = 2 Hz, which gives using equation (TTOl) the 
fractional error of the false alarm estimation 

l 

^FA _ n n ( 1 , n -A 2 __ r .. 







3.2 ^ 






R 


0.32 




0.02 





K-\ - 3 - 2 {r + 0A ) p = 5xl0" y , g = 0.1 

~ 1.05 for ^=100 

~ 1 for R — > oo. 

So the error saturates at 100%, and this saturation is achieved for R > 10. 

Those two examples show that the maximal number of useful time slides and the 
false alarm estimation precision strongly depends on the relative properties of the two 

|j As noted in section [2T4l coincident triggers are denned as less that ±\ At apart. 
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detectors. In particular when there are much more triggers in one detector than in the 
other, the background can be badly estimated and increasing the number of time slides 
does not solve the issue. 

6. Conclusions 

We have studied the statistical error in the background estimation of event-based GW 
data analysis when using the time slide method. Under the assumption of stationary 
noise we analytically computed this error in both the two-detector and three-detector 
case, and found excellent agreement with Monte Carlo simulations. 

The important resulting consequences are: the precision on the background 
estimation saturates as a function of the number of time slides, this saturation is most 
relevant for detectors with a very different trigger rate where the background estimation 
precision can be poor for any number of time slides. 

Let us note that the time slide method can be used in other situations than 
GW data background estimation. For example it can be used to estimate the rate 
of accidental coincidences between a GW channel and an environmental channel in a 
GW interferometer; or in any experiment where coincidences between two (or more) 
trigger generators are looked for. The results of this paper can be straightforwardly 
extended to such an experiment. 

Another limitation, the non stationarity of the data, has not been investigated in 
this paper. Data non stationarity is a well known issue in GW data analysis [12]. In 
the context of the time slides method it raises the question whether the time shifted 
data are still representative of the zero lag data, when large time shifts are used. It 
involves both the problem of the measure of the level of data non stationarity, and the 
estimation of the error it induces on the background estimation. Further work on this 
issue will be the subject of a future paper. 

Appendix A. Two-detector integral 

To compute the integral 

-I J" J 

J xi Jx N Jy 1 Jy N 

-jv ^ Xk < p^ Xi < pMvmv < qMv^q) < v)> ( AA ) 

71"! 7T2 k I 

we put the sums outside the integrals. When k ^ I, the integrals on Xk and x\ are 
independent, and the integration over x±,...,xjf gives a p 2 contribution. Otherwise 
the integration gives a p contribution. Analogously for the y variables we get q 2 or q 
depending on whether TTi(k) ^ 7r 2 (Z) or not. 



7T k 
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Thus we 


get four types of integrals 








integral 


x number of such integrals 


k = Ijir^ 1 


O TTi(k) = I 




x NR 


k ^ l,^ 1 


O TTi(k) = I 


1 2 
JV Vq 


x NR(R - 1) 


k = Ijir^ 1 


O TTi(k) 7^ I 


1 2 
# W 


x JVi2(i2 - 1) 


k ^ l^ 1 


O 7Ti(fc) 7^ I 


1 2 2 


x N[R(R- l)(N-2) + R(N 



(A.2a) 
(A.2b) 
(A.2c) 
(A.2d) 



Here we used that the composition of two circular permutation is a circular permutation, 
and that the only circular permutation with a fixed point is the identity. 
The details of the combinatorics are as follows. 

• k = Ijir^ 1 ° fti(k) = I : There are N different k values. For each of them there is 
only one / that is equal to it. Here it 2 l o tti is a circular permutation with a fixed 
point, so it is the identity. There are R different tti, and for each of them only 
^2 — tti gives the identity. 

• k 7^ I, ix 2 x o iri(k) = I : There are N different k values. For every pair ~K\ 7^ n 2 we 
get 7rf 1 o n 2 (k) 7^ k. And the choice of this pair determines uniquely an / that is 
not equal to k. There are R(R — 1) such pairs. 

• k = l,^ 1 o m(k) 7^ I : There are N different k values. The value of I is 
determined by the equality k — I. And there are R(R — 1) pairs of 7r 1 ,n 2 such 
that Tii 1 ° ^(k) 7^ k. 

• k 7^ I, ix 2 x 7Ti(k) 7^ I : There are N different k values. In the case where 7Ti 7^ n 2 , 
we need that / 7^ k and / 7^ ix 2 X ^i(k), there are N — 2 such /. In the case where 
7Ti = 7T 2 we get k = n^ 1 o 7Ti(k), so there is only one inequality on /, and there are 
N — 1 possible I. 

By summing the 4 terms above and subtracting Mean 2 we obtain 
Var = ^Npq [l + p(R - 1) + q(R - 1) +pq((R - 1)(N - 2) + (N - 1))] - (Npq) 2 



R 

= Npq 
~ Npq 



1 , . pq 



R +p + q 



(p + q) 



R 



2pq 



(A.3) 
(A.4) 

(A.5) 
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Appendix B. Three-detector integral 

We want to compute the integral 



M 9 , , 

R 2 

[XYX'Y' + XZX'Z' + YZY'Z' + AXYZX'Y'Z' + IXYX'Z' + 2XYY'Z' 

+2XZY'Z' - AXYX'Y'Z' - AXZX'Y'Z' - AYZX'Y'Z'} , (B.l) 

where the ' denotes whether the hidden variables are 7i"i, k or 7r 2 , I. 



Similarly to Appendix A we have here eight kind of integrals. 



X 




Y 




z 


number of such integrals 


k = l, 




o 7Ti (fc) = I, 


7T 2 O 7T^ 


\k) = l, 


NR 


k = l, 




o 7Ti(fc) = I, 


7T 2 O 7T^ 


\k)^l, 





k = l, 




o TTi(k) 7^ I, 


7T 2 O 7T^ 


\k) = l, 





k = l, 




O 7Ti(k) 7^ /, 


7T 2 O 7r[" 


\k)^l, 


NR(R - 1) 


k^l, 




o 1Ti(k) = I, 


7T 2 O 7[~ 


\k) = l, 





k^l, 




o 7Ti(/c) = Z, 


7T 2 O 7T^ 




NR(R - 1) 


k^l, 




o 7Ti (ft) 7^ Z, 


7T 2 O 7T^ 


\k) = l, 


NR(R - 1) 


k^l, 




o 7Ti(fe) 7^ Z, 


7T 2 O 7T^ 


\k)?l, 


iVfl - 1)(JV — 3) + (iV- 



In these combinatoric computations we need to assume that all translations are smaller 
than N/4, to ensure that n^ 1 o m o tx 2 x o 7T%(k) = k =>- 7Ti = 7r 2 . This assumption is 
really reasonable, and the result would not be significantly different without it. 
The final result is 



N f 

M 2 = — \ (pq + pr + qr — 2pqr) 
R I 

+ (R— 1) [pq{p + q + pq) + pr(p + r + pr) + qr(q + r + qr) + 6pqr — 4pqr(p + q + r)} 

+ [(R - 1)(N - 3) + (N - 1)] (pq + pr + qr - 2pgr) 2 j, (B.2) 

Appendix C. "OR" case for D detectors 



Using the same heuristic as in section 13.31 and 14.31 we compute the variance of the time 
slide estimation method for D detectors in the "OR" case. This heuristic yielded the 
same results as the exact computation for the 2 and 3 detector case, thus we may expect 
it to stay true in the general case. 

As in equation (fl~4|) . the variance is the sum of the normal Poisson variance 



D j-l 



Var 



Poiss 



N 




(C.l) 



and the variance due to time slides. 
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The estimate of the mean rate is 
Mean = N 



N 



D j-1 

Y Y & + 6p ^ (pj + s Pj 

Li=i i=i 

/ D \ 

Yp< 

VP) / 



D j-1 



j=l i=l j=l 



(C.2) 



(C.3) 



which leads to a variance due to multiple reuse of the data (assuming (5pf ) — — 



N / 



D 



(d \ 



Varsiides/iV =YPi 
3=1 



Y pk 



Yp* 

\&j / Wj ) 
^ D D \ / D 
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This general formula flCTf) is correctly giving back the extra terms in equations and 
( fl4ll for respectively the 2 and 3 detector case. 
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